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1.0 


INTRODUCTION 


1 . 1 Program Background: A Historical Review 

Laser anemometry (or ve Iodine try-as it is now more commonly iden- 
tified in the United States) based on the measurement of the time-of- 
f light of a scattering particle between two laser beams was first dis- 
cussed in 1968 (Ref. 1). A working implementation of the technique was 
made by Dr. R. Schodl in West Germany for DFVLR in or about 1971, and 
his instrument has subsequently been marketed and developed by Polytec 
GmBH. In their systems, the time-of-f light technique is referred to as 
the L2F (Laser-two-focus) principle. 

In 1973 Dr. A. E. Smart, then of Rolls-Royce Ltd in England, 
built a two-spot or transit system based on the same principle, but 
describing the measurement method as laser transit anemometry (LTA), for 
research measurements in aircraft engine high-speed axial compressors. 
The L2F/LTA distinction in names is historical. Insofar as there is a 
difference, it is one of implementation rather than of principle. The 
major dissimilarity is that in L2F the spot source and conjugate 
receiver optical stops are rotated separately by a toothed belt drive, 
whereas the LTA uses an image rotator to spin the image of the test 
space, leaving both sources and conjugate detector stops fixed. The 
first image rotator was a doubled-f olded Abbe prism but was changed to 
one using a mirror Dove design for all subsequent machines because of 
its superior performance with respect to aberrations, transmission 
efficiency, physical and optical length, and reduced weight. Other 
major differences were restricted to the way in which the signal was 
processed. 

After joining Spectron Development Laboratories in California in 
1977, Dr. Smart redesigned his LTA system and pursued further develop- 
ment of the signal retrieval technique, both theoretically and by 
building hardware. Beginning in 1978, the Arnold Engineering Develop- 
ment Center supported the initial phase of development of the "Correlex" 
instrument, designed for ultra-high speed processing of the electrical 


signals from laser transmit anemometers and for other physical and chem- 
ical applications requiring high-speed correlation techniques. In 
succeeding years seven LTA systems, based on the Spectron design of the 
LTA optical head and the Correlex were delivered to NASA, the U.S. Air 
Force and commercial organizations. All the optical heads were 4" lens 
diameter units (Model 104) with a range of optical throws from 367 mm to 
1600 mm (using interchangeable front lens), except for one 6" unit 
(Model 106) which had throws of 1 and 4 meters. 

The initial 5 ns Correlex correlator concept embodied all the 
design features of the later units, but not all were then implemented. 
The first Correlex, appearing as a simple single channel system inter- 
nally based upon the Intel 80/86 microprocessor, was controlled by a 
repackaged Apple computer. Capability for multichannel time gating was 
included from the start, requiring a high-speed external triggering 
system. Two later systems provided this function internally using a 
Spectron proprietary phase address card (PAC), which was itself a real 
time multi-channel sequencing computer. It proved too complex and 
expensive for all but the most extreme channel switching rates, as, for 
example, those required in high-speed multi-stage axial compressors. 
The last machines delivered were controlled by a PDP 11/23 computer, 
which had basic software for controlling the optial head, data acquisi- 
tion, display and storage, and some simple nearly-real-t ime signal 
retrieval processing. 

The Spectron LTA system offered some advantages to complement the 
already extensive fringe anemometry measurement capabilities at NASA and 
Air Force Research Centers. LTA offers high accuracy even from small 
high-velocity particles close to walls, at the expense of measurement 
time. It has good capability for accurate speed and direction measure- 
ment and also for determining the long-term two-dimensional properties 
of statistically stationary turbulence below about 13%. The application 
of the LTA instrument to various fluid dynamic situations has shed much 
light upon the instrument properties and capabilities. These experi- 
ments also demonstrated the need for careful interpretation of the 


derived measurements and an enhanced understanding of what is now per- 
ceived to be a highly complex technology. The detailed study and com- 
parison of data obtained from fringe aneraometry and LTA techniques has 
also given insight into both measurement technologies, as well as the 
fluid dynamic properties of the flows under examination. This experi- 
ence has facilitated development of the LTA data processing and informa- 
tion retrieval software now available at NASA and Air Force Research 
Centers, and has established greater confidence in the results obtain- 
able from LTA systems. 

Details of the investigations performed by Spectron for AEDC, 
which provided much of the background to the developments described in 
this report, can be found in Ref. 2. Further information relating to 
the Correlex principles, design and performance is available in Refs. 3 
and 4. 


1.2 Program Goals 

The primary program goal has been to provide working code which 
will enable the LTA systems to approach the performance in real applica- 
tions that theoretical analyses suggest should be possible. To meet 
this objective a number of activities have been undertaken. These 
included the design of three major software packages; for signal acqui- 
sition and system control (ACQUIRE), for data retrieval, processing and 
archiving (PROCESS), and a data simulation program to permit verifica- 
tion of the processing package (SIMULTAC). Other subsidiary tasks 
included coding, testing, debugging, verification, and documentation of 
the above software packages, and the provision of a number of supporting 
packages and auxiliary routines. 

1.3 Program Summary 

Earlier SPECTRON work performed for AEDC between April 1st, 1981 
and December 31st, 1983 was reported on in Ref. 2. The experimental 

equipment, consisting of a f ully-operational SPECTRON LTA System with 
the Correlex as supporting signal-processor, was described in four 
volumes, supplied with the instrument: 


Vol. 1 Optical Head/SDL 81-53003 
Vol. 2 Correlex Operator’s Manual/SDL 82-53010D5 
Vol. 3 MSCLTA II User’s Manual, Version 1.1 
Vol. 4 MSCLTA II User’s Manual, Version 2.0 Rev. 


- March 1981 

- October 1982 

- February 1983 

- September 1986 


The present program was initiated to produce software capable of 
exploiting the full potential of the SPECTRON LTA System. A major 
activity was the establishment of a modular architecture within which 
efficient software development could take place. The importance and 
validity of this effort are clearly demonstrated by the relatively minor 
changes that have been found necessary during the course of the program. 
VAX PASCAL was chosen as the programming language, and the specified 
target hardware was the VAX 11/730 with VMS operating system. The soft- 
ware structure having been established, a large amount of work was 
devoted to error trapping, documenting, and rendering robust the code 
for ACQUIRE and PROCESS. Once this development vehicle was available 
algorithm development proceeded in the confident expectation that 
programs would not crash because of unexpected errors in the house- 
keeping and organizational code. PROCESS, as it was developed to 
include more of the functional modules, began to reveal problems in 
running the initial simulation program DAS. As a diagnostic tool, DAS 
was found to have insufficient power. A new program, SIMULTAC, was con- 
ceived and a preliminary version implemented. 

The status of the programs delivered at the end of the contract 
period is as follows: 

ACQUIRE is conplete, robust and fully operational. 

SIMULTAC is as complete as has been found possible with the resources 
available. There are always, of course, elaborations and enhancements 
that can be added to simulation, test and diagnostic tools. 

PROCESS is complete, robust and fully operational. The embedded 
algorithms have not been evaluated over a sufficiently wide range to be 


sure that they are optimized for their current purpose. More work on 
the structure of these algorithms and their implementation should result 
in more powerful and efficient operation, and might make better results 
available from poorer data. A framework has been constructed within 
which such improvements may evolve. Specific recommendations can be 
found in Section 4. 

The key objective of the effort was the development in an optimum 
form of algorithms and related computer codes, together with appropriate 
documentation, for extracting from the LTA system correlograms obtained 
in a general three-dimensional flow field the two components u and v of 
the mean velocity, the root mean square fluctuations u 1 and v f in these 
components, and the first-order velocity shear stress u'v f « 

A pair of matched discriminators, designated Channel 8 and having 
1.6 microsecond detection filters and, 6.4 microsecond center-estimation 
filters, were to be delivered to Langley for installation and eval- 
uation. A second pair, designated Channel 7 and having 800 nanosecond 

detection filters and 3.2 microsecond center-estimation filters, were 
also constructed and delivered to Langley. Hie characteristics of the 
discriminator modules available in the Langley instrument are discussed 
in Section 2.5. 


2.0 


DEVELOPMENT TASKS 


2.1 Data Acquisition and System Control Software - ACQUIRE 

In the preliminary versions the entire program for operating the 
LTA system existed as a single suite of modular code. The decision was 
taken at an early stage to separate ACQUIRE and PROCESS into two auton- 
omous programs which could be run concurrently within a multi-tasking 
environment. This option enabled ACQUIRE to be run in real time with no 
impairment of data acquisition rate, while simultaneously making avail- 
able selected unprocessed correlograms at the CORRELEX display. System 
parameters are continuously displayed at the command terminal. 


2.2 Simulation Software - SIMULTAC 

A major effort was devoted to the further development and refine- 
ment of the LTA simulation program SIMULTAC. The program is optimized 
for speed, using deterministic analytically-based computations to gener- 
ate the basic correlograms. The effects arising from a finite popula- 
tion of light-scattering particles are simulated by perturbing the 
correlogram, using a Poisson random number generator. Background noise, 
assumed to be additive, is generated in the form of appropriately scaled 
Poisson distributed random numbers. 

The basic physical assumptions underlying both the simulation and 
data-processing software are that: 

(i) the particle’s velocity is constant between the light spots; 

(ii) the particle concentration is constant throughout the fluid in 
space and time; 

(iii) there is negligible relative motion between particle and fluid; 

(iv) the flow is statistically stationary; 

(v) the optical geometry and light intensity distributions are 
unaffected by the flow properties. 


2.2*1 Theoretical Background to SIMULTAC 

The geometrical model on which SIMULTAC is based (see Figure 1) 
assumes two identical light distributions, whose central maxima 0^ and 
are separated by a distance s. It is further assumed that the two- 
dimensional light distribution remains constant in the axial direction, 
at least over the region of experimental interest* O^X is the reference 
axis, in the plane normal to the axis of propagation, from which the 
orientation a of the line joining the spots is measured. 

The circles represent spatial contours within which the scattered 
light level from a given particle exceeds the threshold of detectabil- 
ity. At present the threshold circle is assumed to have some fixed rad- 
ius, q; this is equivalent to making the additional assumptions that all 
the scattering particles are of the same size and that all particles 
which pass within a distance q of 0^ or 0^ are detected, whatever the 
speed. (A modification which would take into account the fact that the 
light distribution is in reality Gaussian and that faster particles 
scatter less light is discussed later.) u is the velocity component in 
the direction O^X; v is the orthogonal component in the plane of the 
diagram. 6 and 0 are the angles made by a general particle trajectory 
with O^X and 0^0^ respectively. Hence: 

9 = a + 0 


For a particle transit to contribute to the correlogram the cor- 
responding trajectory must intersect both threshold circles. In addi- 
tion, the transit time (time between points of closest approach to the 
propagation axes) must be less than the maximum displayed correlogram 
lag. 

We can now calculate the bounds imposed by the threshold con- 
straints on the permissible values of 9. 

For a given 9 there is a band of width w, say, disposed symmet- 
rically about the line O^C^, which defines those trajectories which make 
non-zero contributions to the correlogram. 
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Figure 1. Transit geometry 






The extreme trajectories in this band are tangent to one or other 
of the threshold circles, and from Figure 1, 


w = 2q - s |sin B| 

or 

w = 2q - s |sin (0 - a) | . (2.1) 


Note that w/2q may be interpreted as a probability that an event at the 
first spot is followed by an event at the second. The maximum permis- 
sible absolute value of 3 is determined by the direction of the particle 
trajectory for which w becomes zero. Thus, 

sin 3 = 2q/s , 

max 

and hence 


ol - sin“l (2q/s) < 0 < a + sin"l (2q/s) . (2.2) 

This then determines the range of 0 within which a particle trajectory 
may make a contribution to the correlogram. 

The next step is to define the probability density distribution 
for the velocity components u and v. This function is assumed to be a 
joint Gaussian random process of the form 

< \ = 1 r 1 r (u ~ u) 2 

^uv u>v 2ttct a /l-p z ex P L 2(l-p : 0 * 

x y x 

(2.3) 

2p(u ~ u)(v - v) . (v - v) 2 i i 

+ rsru > 


where TT and "v are the mean values of the velocity components u and v 
respectively, and cr^ and q ^ are the corresponding standard deviations, 
p is the correlation coefficient, which determines the interdependence 
of u and v, and is associated with shear stress through the relationship 


pa a 
x y 


(2.4) 


where u 1 and v 1 denote the fluctuating parts of u and v. 


In the early versions of SIMULTAC this continuous probability 
density function was reduced to a two-dimensional histogram of values 

corresponding to discrete cells (i, j) in (u, v) - space. For the 

(i,j)th cell 

9 = tan * (Vj/u^) 

and only cells for which this value of 9 fell within the range given 
above contributed to the correlogram. 

A major improvement in speed and efficiency was achieved by 
carrying out the computation of the correlogram directly in (x, 9) - 

space. The relationship between the probability distributions in the 
two spaces is 

P T0 (T ’ Q) = ~o P uv (u> v) (2 ‘ 5) 

(A table of transformations between various probability distribution 
spaces is given in Section 2.3.4.) 

Then the content of the ith bin of the correlogram at angle 

setting otj is computed as 


3 1 T. - 0. 


Ti + 0.5 
5 


aj + A0 i 

J w(0)p uv (u - ^ cos0, v = “ sin0)d9dx 


( 2 . 6 ) 


where 

A9 * sin 1 (2q/s ) , 

w( 0) = 2q - s | sin ( 9 - a) | > 

and the assumption is made that contributions to the ith bin arise from 
events whose associated delay times are symmetrically disposed over one 
sample time around . 



Note that an extra factor x- 1 appears in the integrand of 
Equation (2.6), which represents the effect of velocity - biasing. This 
formula for in which the width w of the trajectory acceptance 

band appears as a weighting function, can be established by the 
following argument. 


2.2.2 Dependence of the Correlation Function on the Flow-Angle 

Weighting Function and Velocity 

We assume that the flow is two-dimensional in the plane normal to 
the beam axes and that the detector fields of view cover an axial 
distance £• 

Suppose the flow direction is in the range (6, 0 + 60) and the 
flow velocity is in the range (V, V + 6V). Then transit times are in 
the range (x, x + 6x), where x = s/V and s is the spot separation. The 
probability of occurrence of these conditions is 

>xq (t, 0) 6x 60 


, .v 5 ~ ( S COS 0 s 3ino\ 

(X, 0 = 1 p lU = , V = 

X° UV v X X J 


s sin9^ 


X0 


ective cross-sectional area presented to the flow is w£ where, as 


w - 2q - s | sin ( 0 - a) | 

(.s the threshold radius. 

t If the particle concentration is c, the number of particles 
jag this area per unit time is, for the given flow conditions, 

s 

j cw£V = cw£ - 
x 
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2.2.3 Program Description 

Initial input parameters include the mean flow velocity and 
direction, the axial and cross -turbulence intensities (standard devia- 
tions as percentages of the mean values), and the correlation coeffi- 
cient, p. The number of standard deviations can also be specified at 
which computation of the correlation function is ended, and all 
remaining channels are set to zero. The purpose of this feature is to 
economize on computational time at very low turbulence levels. The LTA 
system parameters which must be set include spot separation and spot 
diameter, threshold radius, spot orientation (relative to mean flow), 
angular increment, number of stores, delay time and acquisition time per 
correlogram ( f stop time 1 ). 




The simulation then begins by computing a set of correlation 
functions according to Equation 2.6, which includes the velocity-biasing 
factor. In the interests of accuracy, the double integral is evaluated 
for each correlogram bin by Gaussian quadrature, using repeated one- 
dimensional integration (see, for example, Ref. 5). There are 10 nodes 
in each integration unless o ^ or is less than 1%, in which case 
20 nodes are used to ensure adequate sampling of the probability distri- 
bution function. For some test purposes, a zero threshold may be 
required. Then w( 0) is replaced by a delta-function and the double 
integral becomes an integral over t alone. The necessary modification 
to the integration routine is automatically performed in SIMULTAC when a 
threshold value of zero is chosen. 


The operator is then asked whether the correlograms are to be 

scaled to conform with a specific requirement for the highest channel 

content, denoted by h^, of the central correlogram. If not, the program 

stops, the output at this point consisting of a stored set of noiseless 

correlograms. The channel contents of each correlogram are non-integral 

real numbers. If the correlograms are to be scaled, a value for h is 

P 

requested, together with a value for the ratio of mean background noise 

level h to h . This background noise is assumed to be flat in the 
n p 

mean, and represents mainly the effects of parasitic radiation in an 
experiment. (At particle concentrations normally encountered, it can be 
shown that noise-particle and interparticle correlation events are rela- 
tively rare. See also the comments in Section 2.3.3 concerning non-flat 
backgrounds • ) 


The next stage of the program consists of scaling the correlo- 
grams so that the highest peak takes the value h - h ; at this point, 

P ^ 

the correlograms contain contributions from particle events only. The 
operator is then asked if f finite seeding’ is to be simulated; i.e., if 
statistical effects equivalent to those resulting from a limited number 
of scattering centers are to be included. If not, and if a zero value 
has been entered for h^/h^, execution of the program stops. If finite 
seeding is required, the correlograms are randomized by means of a 


Poisson number generator. The assumption underlying this step is that 
the particles are purely randomly distributed in the fluid. Then the 
number of contributions during the experiment from particles within the 
speed range (V, V + 5V) will be a Poisson-distributed random variable. 
For each correlogram bin, the mean of the Poisson process is set equal 
to the bin content. It should be noted that, since a Poisson process 
involves integer random numbers, a rounding up or down of the bin 
contents is necessary. 

The total number of particles contributing to the entire correl- 
ogram set can now be calculated. Since each correlated particle contri- 
butes a f l f to the correlator output, this number is simply the sum of 
the correlation bin contents over the complete set of correlograras . 

If h^/h^ is zero, execution of the program stops at this point, 

the output consisting of a set of correlograms in which the only noise 

source is the finite number of contributing particles. If h /h is non- 

^ P 

zero, a randomized noise background is added to each bin. This is taken 
to be a Poisson process with mean h^. The correlograms now represent 
the result of an experiment involving a finite number of particles and 
an additive noise background with constant mean. The operator is then 
asked if jitter noise is to be included in the experiment. This is 
intended to simulate errors in transit time estimates arising from, for 
example, low light levels or electron transit time spread or variable 
gain in the photomultiplier tube, and is only significant at very low 
turbulence levels. The errors are modeled as a Gaussian random process 
and are under most experimental circumstances insignificant compared 
with other noise sources. For this reason, and by mutual agreement, 
development of this part of the program was given low priority and, 
although the necessary hooks are in place, coding of the procedure 
JITTER is largely incomplete. The intended design consisted of a 
Gaussian randomization of the delay time associated with each count in 
the correlogram; the resulting random jitter will occasionally move 
counts into adjacent bins. The standard deviation for the Gaussian 
process will depend on the discriminators in use (see Section 2.5). 
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Finally, differing correlogram seeding rates can be simulated. 
If the ’seeding change 1 option is selected, each correlogram is scaled 
by a different random number generated from a Gaussian distribution with 
a mean of unity and a standard deviation set by the operator. In order 
to restrict fluctuations to a realistic range, only random numbers with- 
in one standard deviation of the mean are accepted by the program. This 
simulation of fluctuating seeding rates is termed ’denormalization’ . As 
described below, PROCESS includes a compensating 'renormalization' 
procedure. 


2.2.4 Comment on Threshold Modeling 

The assumption of a fixed threshold for particle detection in 
Section 2.1.1 is clearly unsatisfactory. A rather more realistic model 
can be constructed in the following way. 

We retain the assumption that all the particles are of the same 
size, but now assume that the light distribution in each spot is 
Gaussian. Since the total quantity of scattered light is inversely pro- 
portional to the particle speed V, we shall assume, as a first-order 
approximation, that the peak of the output from the Gaussian discrimina- 
tion filter is also inversely proportional to V. A particle transit is 
detected if this peak exceeds the discriminator threshold setting. To 
take account of the Gaussian light distribution, a factor 
exp f- 2r 2 /r 2 ], where r is the beam radius and r is the distance to 
the beam axis at the point of closest approach, must be included in the 
calculation. Then we can write for the peak height I of the filter 
output 

O 

where k is some constant of proportionality. 

Suppose the threshold setting is such that, at the mean speed V , 
a particle is detected when r^ = q. Then the criterion for a detection 
to occur at any other speed is 
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2.3 LTA Data Reduction Software - PROCESS 

The following discussion is concerned with the mathematical basis 
for the sequence of data reduction operations which appear within 
PROCESS and whose ultimate objective is the extraction of the five flow 
parameters, consisting of mean velocity components u and v, their r.m.s. 
fluctuations u f and v f and the first order shear stress u f v f . The main 
stages in the reduction of the data are: 

1. Renormalization of correlogram set 

2. Correlogram filtering 

3. Correlogram background removal 

4. Transformation of the data to (V,9) space 

5. Extraction of flow parameters 

2.3.1 Renormalization 


Variations in the rate of acquisition of data from one angle 
setting to another can result in a partially random scaling within a 
given set of experimental correlograms . Simulation of this feature of 
the physical LTA environment, referred to in the program comments as 
denormalization, is available as an option in SIMULTAC. 

To compensate for these data rate fluctuations, processing of the 
data currently begins with a renormalization stage. The peak correlogram 
value y^ for each angle setting 0^ is located and a best fit found to the 
Gaussian formula 
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where c, 0 and n are unknown constants. Each correlogram is then 

o 

rescaled by the new peak value. Note that the assumption that the 
points lie on a Gaussian curve may not be strictly correct. Numerical 
experiments with simulated data, however, suggest that the approximation 
is a very good one. 


2.3.2 Correlogram Filtering 

A self-adaptive digital filter is used to smooth the raw data. A 
triangular (Parzen) weighting function is used, which has better peak- 
preserving properties than a rectangular shape and, hence, provides 
better performance at low turbulence levels. The background level for 
each correlogram is first estimated by averaging over the data between 
the third store and the store nearest to 0.6 times the peak store 
number. After subtracting this average level the remaining data are 
summed, giving an estimate of the total signal content, Q say, of the 
correlogram. If H is the peak store content, the half -width of the 
filter is defined as 

HW - ROUND (B > 

where the operation ROUND makes HW an integer. 6 is an empirical con- 
stant which provides control over the broadening induced by the filter. 
Any filtering, of course, involves a compromise between noise reduction 
and loss of signal detail. Further investigation of the relationship 
between turbulence level and the optimum value for 6 is clearly desirable. 


Having determined HW, the filter is now applied to the original 
correlogram data. Each data point Y(k) is replaced by a triangularly 
weighted average over the set of data points from locations k - HW to 
k + HW according to the formula 


Y(k) * 


) Y(j). 



= k + HW 

I 


= k - HW 


(1 



The filtering computation is carried out in the procedure FILTER 
in the ALGORITHM module. 

2.3.3 Correlogram Background Removal 

During the course of this work much discussion has been devoted 
to the problems of dealing with background noise in the correlation 
data. Unfortunately, some characteristics found in experimental data 
are not fully understood; for example, background levels may differ 
significantly at small and large delay times, and sloping backgrounds 
have been observed on a number of occasions. The basic assumptions are 
made in PROCESS that the background can, for most purposes, be reason- 
ably well approximated by a horizontal correlogram floor, and that the 
signal content does not occupy more than the central 50% of the correl- 
ogram. Accordingly, the means of the contents of two groups of 

channels, occupying the initial and final 20% of the correlation 
function, are calculated and compared. The greater is taken as the 
effective overall noise level. The channel numbers at which the central 
correlogram peak decays to the noise level are also determined at this 
stage. After subtraction of the background, the correlogram is 

truncated by setting all channel contents to zero beyond these points. 
The computations are performed in the ALGORITHM module in procedure 
BACKGROUND. 

If, in a particular experiment, the correlogram background is not 
flat, the effect of these numerical procedures will be to over- 
corapensate for the noise and reduce the signal contribution. For 
unimodal velocity distributions, mean values will probably not be 
greatly affected, but higher moments will certainly become unreliable. 
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2.3.4 Transformation of the Data to (V, 0) Space 

To clarify the discussion, Table 2.1 is included here to show the 
relationships between the various spaces in which the velocity probabil- 
ity density function can be expressed. These relationships are governed 
by the modulus of the Jacobian of the appropriate transformation, |j|, 

TABLE 2.1 


FUNCTION SPACE 

INPUT 

OUTPUT 

(u,v) 

(V , 0) 

(u , v) 

(t, 0) 

(t,0) 

(V,0) 

(t, 0) 

(u,v) 


PDF 

RELATIONSHIPS 


S 2 / T 3 

t 2 /s 
t 3 /s 2 


P V 0 Vp uv 


P T 0 3 P uv 


P V 0 S P T 0 

T 3 

P uv 2 
s z 


The variables in the table have their conventional meanings and 
are related by the formulas 

u = V cos 0, 

V = /u 2 + v 2 , 
t = s/V. 

It is appropriate at this point to correct an error appearing in 
Ref. 2; on Page 33, the factor t 4 should be replaced by t 3 in Equations 
(4.7) and (4.8). 


v = V sin0, 

0 = arctan (v/u), 
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All the flow information is, of course, available in any of the 
spaces of Table 2.1. The choice of space in which to carry out the 
data-processing will be governed in practice by a number of considera- 
tions, including noise characteristics, the significance of nonuniform 
spacing of samples, interpolation or truncation problems and computa- 
tional burden. It is not possible, on the basis of present experience, 
to make very general recommendations, but a scheme which operates almost 
entirely in (u,v) space and which is based on weaker assumptions than 
the present procedures is described in Section 4.1. 

Velocity-biasing should be considered in the treatment of 
experimental data. As demonstrated above (see Section 2.2.2) for the 
simple fixed-threshold detection model, the unbiased distribution func- 
tion can be recovered by including an additional division by velocity, 
or, equivalently, multiplication by t/s . (See also Ref. 6) For 
example, the transformation from ( t, 9) space to (V,0) space is given by 


P V 9 


(V, 9) 



( 2 . 8 ) 


If p^ represents the biased distribution function obtained from experi- 
ment, we then have 


P V9 (V ’ 9) = Z* 


T 0 


( T, 9 ) 


(2.9) 


In a Correlex-LTA system, each detected transit is represented by 
a single pulse at the appropriate delay time. Hence, in the absence of 
noise and given adequate sampling of the flow, the final correlogram 
should have the same form as the biased distribution function expressed 
in (t, 9) coordinates; i.e., it should be proportional to p' . Thus, a 

T 0 

section through the p distribution can be recovered directly from the 

v 0 

correlogram by the use of Equation (2.9). 



In the current version of PROCESS, the correlograms are transformed 
from (t> 9) space to (V,0) space before the required flow parameters are 
extracted. The necessary transformation of the "tr-coordi nates into V-coordinates 
is accomplished using the formula 


i v = TRUNC [ 


16641 
i + 1 


1 ] 


T 

where i-f is the correlogram channel number and iy the corresponding channel 
number in V-space. Note that 16641 = (256 x 65) + 1. The '1' is added to 
ensure that a significant error is not introduced when the divisor is a factor 
of 16640. Only channel numbers 64-225 are converted in this way (into channel 
numbers 255-64 in V-space), the assumption being that there are no signal 
contributions to the correlogram for i ^ = 0-63. The data in V-space, now non- 
uni formly spaced, are also taken to be zero for iy = 0-63. 


2.3.5 Extraction of Flow Parameters 


The correlation function in (^T.O) space has now been transformed into an 
unnormalized section through the velocity probability density function in (V,0) 
space. The succeeding stage of data-processi ng consists of extracting from each 
of the transformed correl ograms , at angle settings 8 j , the peak height Hj , peak 
position P j , the sum of the channel contents Qj and the first and second moments 
of velocity Vj and Tj , respectively. The best parabolic fit is then found to 
the 5 highest values of the set ^HjJ, using the equation 

H = 3l e 2 + bjO + c l (2.10) 

At the highest point (H 0 ,9 0 ) of this curve dH/d0 = 0, giving as the most prob- 
able flow angle 
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9 = 


- y2aj 


Then H - c, - 
o 1 


1L 

4a, 


This calculation is based on the assumption that the velocity 
distribution can be well approximated near the peak by a parabola. This 
is not unreasonable for a joint Gaussian function. We shall also assume 
that B q is a good approximation to the mean flow angle. 

Having found 0 q , the nearest five angle settings can be identi- 
fied. A second parabolic fit is now performed on the five values of 
associated with these angles, using 


V = a 2 e 2 + b 2 e + c 2 . 


( 2 . 11 ) 


Having thus determined a 2> t> 2 and c 2> the mean speed is taken to be 


We now have our estimates of mean flow speed and direction. 
Estimates of and P are obtained by locating pairs of points in 

each transformed correlogram (still in (V,8) space) which lie on the 
contour of the velocity distribution function at some given level. In 


the current version of PROCESS this level is l/\Te* times the highest 
point H q of the distribution. 

This choice was originally made because the slope of the proba- 
bility density function is a maximum at this plane, and the sensitivity 
of the data points to noise is minimized. Selecting a lower level, how- 
ever, has the advantage of generating more intersections, and the choice 
is not necessarily optimal. 
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The basic assumption here is that each correlogram in (V, 0) space 
can be described by a Gaussian. Although not strictly true if p is 
significantly non-zero, the errors are probably very small for contours 
at the l//e level. Hence, we write for the jth transformed correlogram 


, C v “ V i ) 2 , (2.12) 

Y = H exp { 2T ^ } 

j 


The l//e plane is at an absolute height of H // e . Setting Y 

o 

equal to this value, we obtain for the pair of points on the contour 


v e j = V j - T j I ” 2 l0 S (H o /H j 

Pairs of points can be obtained in this way from all correlograms 

for which H. > H //e. The coordinates (V .,0.) are finally transformed, 
jo ej j 

also in procedure KARY, to (u,v) space by means of the relations 

u . = V . cos (0. - 0 ) , 
ej ej jo 

(2.13) 

v . = V . sin (0. - 0 ) . 
ej ej jo 


Referring the angular settings to 0 q aligns the v-axis with the mean 

flow. The derived values for a and a will now represent the required 

x y 

axial- and cross-turbulence magnitudes respectively. 

It should be noted here that there is an error in the procedure 
described above. The pairs of points at the l//e level should be found 
in (u,v) space, not (V, 0) space, since the transformation from one space 
to the other alters the shape of the probability distribution function. 
The appropriate Jacobian is given by Table 2.1. 


The parameters a^, and p are obtained by the following 
argument. 


Equation (2.3) represents the model for the normalized velocity 
probability density function. An arbitrary scaling factor is introduced 
in the case of the function, p 1 say, describing the experimental data: 
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. , r 1 r ( u " u) z „ 

P U v - h o ex P L " 2<lV) i — r* 2 


(u - u) 2 n (u ~ U ) ( V ~ v ) 

p 


a o 

x x y 




(2.14) 


At the l//e contour level 



88 h //e . 
o 


Hence the equation for this contour, an ellipse, is 


(u 



2p 


(u - u)(v - v) 


a a 
x y 



Differentiating, we find 


(2.15) 


dv 

du 


(u - u) - 


pa a 
xy 


pa a 
x y 


(u - u) - 


(v - v) 

1 (v - v) 


dv 

du 


00 when 


pa a 
x y 


(u - u) = a 2 (v - v) 
x 


Substituting in Equation (2.15) we obtain 


u 


u + 


a 


x 


When 


dv 

— - 0 we obtain 
du 


v = v + a 

y 


(2.16) 


(2.17) 


Hence a and a can be found from the points on the ellipse at which the 
x y 

tangents are parallel to the coordinate axes. 
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When v = v. Equation (2.15) gives 


p = [ 1 - 1 1/2 (2.18) 

o 2 

X 

Hence p can be found from either of the points of intersection of the 
ellipse with the u-axis. As described above (see Equation (2.13)), the 
coordinate axes are chosen so that v = 0. 

The final steps in the present procedure for computing o , a and 

x y 

p then consist of a least squares fit of an ellipse to the set of points 
given by Equation (2.13), finding the points at which the tangents to 
this ellipse are parallel to the axes, finding a point on the u-axis at 
which v = 0, and finally applying Equations (2.16), (2.17) and (2.18). 


If a contour is chosen at some height other than the e"l/2 level, 
say at e"b 9 it can easily be shown that Equations (2.16) to (2.18) 
become 


(u - u) 2 a 
(v - v) 2 - 
p - [ i - 


1 2bo 2 

X 

when 

dv = 
du 

(M 

O 

CM 

when 

r - 0 

du 

(u - u) 2 - 1 1/2 
2b o 2 

X 

when 

< 

ii 

< i 

• 

the comment may 

be 

added that 


a , a and p can be 
x’ y 


derived directly from the coefficients found in carrying out the ellip- 
tical fit to the data. The equation on which this procedure is based is 


CjU 2 + C^v 2 + uv + c^u + c^v = 1 


By expanding Equation (2.15) and comparing coefficients, we find 


o=2 

x 


c 2 (l 


+ CjU 2 ) 

T 


4 c 1 c 2 


- c. 


1/2 


o=2 

y 


c i (1 


+ CjU 2 ) 

T* 


4c 1 c 2 ‘ c 3 


1/2 


P 


C3/2 ^ c i c 2^ 


1/2 


Similar formulas can be found for contours at other than the l//e" level. 

2.3.6 Multi-station Provisions in the MSCLTA II Software 

The original design of the MSCLTA II software included provisions 
for the control of a Correlex system containing a PAC card for multi- 
station data acquisition. To verify that the multi-station provisions 
did indeed work, a simulated test flow was set up at SPECTRON, where 
multi-station acquisition was successfully performed. Communication was 
established and data sets transferred from the Correlex to the VAX and 
stored. 

Following the first software delivery to the NASA Lewis facility, 
the Correlex equipment was set up. NASA Lewis provided a test flow to 
demonstrate the functionality of both single and multi-station acquisi- 
tion modes. Once again, data were successfully acquired. However, some 
communication problems were encountered with the Lewis arrangement that 
were not detected with the SPECTRON configuration, primarily in the 
setting of certain PAC card front panel parameters. 

While the PROCESS portion of the MSCLTA II software was being 
developed, the simulation program was modified to produce correlogram 


1 
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data sets imitating a multi-station acquisition mode. This permitted 
testing of multi-station features in the PROCESS program.lt is believed 
that the multi-station mode functions as intended, and demonstrated at 
the Lewis facility, although the analysis of data using PROCESS in the 
multi-station mode has not been performed as extensively as for the 
single-station mode. 

2.4 Software Testing 

A restricted set of software tests has been conducted. A strat- 
egy was developed with the goal of achieving a working SIMULTAC/ PROCESS 
combination which met certain performance criteria. SIMULTAC and 
PROCESS were demonstrated initially on ideal data and subsequently with 
a limited number of contributing particles and with background noise. 
Performance criteria included the following: 

(i) Both computers involved - SPECTRON's VAX 11/730 and IBM PC/AT - 
should yield the same results with common software. 

(ii) Output estimates for V , 9, a , a and p should improve in 

x y 

accuracy with increasing numbers of particles, reducing turbu- 

lence and reducing noise levels. 

(iii) Repeated trials with noise data should reveal no biases. 

The test parameters used for (i) and (ii) are set out in 

Table 2.2. Some statistical tests were also carried out for the two 

turbulence levels and for h = 500, h /h “ 0.1 . Ihe results of all 

P n P 

these trials are included as Tables 2.3, 2.4 and 2.5. It can be seen 
that, for the IBM and VAX 11/730 computers, the performance criteria 
have been generally well met. 
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Table 2.2. SIMULTAC /PROCESS Test Cases 


Spot Separation 
Spot Diameter 
Threshold Radius 

852 micron 
40 micron 
0, 1, 15 micron 

Delay time 
Number of Stores 

200 ns 
256 

Mean Velocity 
Mean Flow Orientation 
Axial Turbulence 
Cross Turbulence 

30 ms- 1 
0 degrees 

1,5 % 

Correlation Coefficient 

0.4 

Angular Steps 
Angular Increment 
1.1° at 5% turbulence 

9 

0.15° at 1% turbulence 

Height of Correlogram Peak, h 

Fractional Background Level, h^/h 

n p 

500, 100 
0, 0.1, 0.5 
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TABLE 2.5 

STATISTICAL TESTS (VAX COMPUTER) 


Input Parameters: 


Mean Velocity = 30 ms- 1 

Threshold Diameter = 30 ym 

Greatest Channel Content (h ) = 500 

Fractional Background Level ^ = 0.1 

Output Estimates for cr = a„ = 5% 


Run No. 

Velocity Estimate 
(ms- 1 ) 

Flow 

Direction 

a x 

a 

y 

P 

STAT005/1 

30.076 

0.06 

4.98 

5.33 

0.37 

2 

30.017 

-0.02 

5.18 

5.07 

0.35 

3 

30.067 

0.01 

5.07 

5.52 

0.39 

4 

30.064 

0.04 

5.01 

5.23 

0.37 

5 

30.056 

-0.02 

4.92 

5.13 

0.38 

6 

30.055 

-0.03 

4.95 

5.45 

0.38 

7 

30.076 

0.01 

5.01 

5.27 

0.38 

8 

30.093 

0.04 

4.99 

4.91 

0.37 

9 

30.072 

0.11 

5.05 

5.44 

0.39 

10 

30.052 

-0.04 

4.99 

5.48 

0.41 


Output Estimates for <j > = 


Run No. 

Velocity Estimate 
(ms- 1 ) 

Flow 

Direction 

a x 

a 

y 

P 

STAT001/1 

30.001 

-0.12 

1.13 

1.74 

0.25 

2 

29.999 

-0.12 

1.07 

1.98 

0.17 

3 

29.993 

-0.11 

1.15 

2.00 

0.33 

4 

29.999 

0.06 

1.07 

1.65 

0.17 

5 

30.003 

-0.05 

1. 19 

3.05 

0.46 

6 

30.005 

-0.03 

1.13 

1.79 

0.21 

7 

29.991 

-0.04 

1.13 

1.86 

0.27 

8 

30.006 

-0.03 

1.13 

1.68 

0.20 

9 

30.011 

0.05 

1.20 

1.39 

0.14 

10 

29.995 

-0.06 

1.11 

1.95 

0.22 























2.5 


Pjgcripina.tpcs 


As part of the third-year contract, additional discriminator 
hardware was provided to extend the LTA equipment capabilities to 
include lower speed flows* Each of the six discriminator modules on a 
discriminator card, of which there are two in the optical head, is 

intended for a specific speed range, with a factor of two between 
ranges. Each discriminator is designed to be compatible with flow 

speeds from one-half to twice the center value, although the performance 
is acceptable over a range several times greater than this. The six 
pairs of discriminator filter modules originally supplied to NASA 
Langley are described in Section VI of SDL 81-53003 (Optical Head), 
March 1981, which was provided with the original LTA equipment. The 
ranges originally covered were identified as 1 through 6; the ranges of 
the two new filter modules supplied under the present contract are 

designated 7 and 8. The characteristics of the detection filters and 

center estimation filters are shown in Table 2.6. 

Table 2.6. 


Slot 

Range 

No. 

Detection Filter 
Width 

Center Estimation Filter 
Width 

A 

i 


15 

ns 

50 

ns 

B 

2 


25 

ns 

100 

ns 

C 

3 


50 

ns 

200 

ns 

D 

4 


100 

ns 

400 

ns 

E 

5 


200 

ns 

800 

ns 

F 

6 


400 

ns 

1.6 

Vis 

- 

7 


800 

ns 

3.2 

ps 

- 

8 


1.6 

VS 

6.4 

Vis 


Figures 

2-6 

which relate 

to the 

Range 7 discriminators, are 


included to illustrate typical performance features. 
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Figure 3. Impulse approximated by a 5 V square wave 100 ns wide. This 
plot shows the response of the detection filter of Board // 2. 




Figure 4. Impulse approximated by a 5 V square wave 100 ns wide. This plot 
shows the response of the center estimation filter of Board #1. 





Figure 5. Impulse approximated by a 5 V square wave 100 ns wide. This plot 
shows the response of the center estimation filter of Board if 2. 




Figure 6. Impulse approximated by a 5 V square wave 100 ns wide. Total circui 
delay from impulse to output is 2400 ns. This plot shows the output 
delay between Board #1 and Board #2 to be approximately 10 ns. 





A brief description of the operation of the discriminators was 
given in the Manual referred to above. A more detailed discussion can 
be found in Ref. 7. 



3.0 


COMMENTS AND CONCLUSIONS 


3.1 Software Structure 

At the beginning of this program it was realized that success or 
failure in creating a usable set of software depended on a clear under- 
standing of its intended function. Indeed, the process of ’design’ has 
been described as ’understanding the requirements.’ The value of the 
detailed discussions of the framework of the software structure, which 
occupied so much of the initial effort, has been demonstrated by the 
fact that the architecture has subsequently remained essentially 
unchanged. Several modules have undergone some alterations, and cos- 
metic work is still needed, but the skeleton is fixed and has proved 
appropriate for the task. 

The separation of the data acquisition program, ACQUIRE, from the 
information retrieval program, PROCESS, was an early and significant 
development. Its value is obvious in retrospect. The first reportable 
activity consisted of generating a user-acceptable interface; as a 
result, many of the software operations, while accessible to the dedi- 
cated analyst, are transparent to the user in routine applications. 
This early planning made the final programs easy to maintain, and the 
generation of the necessary documentation a straightforward task. 


It should be noted that, in its present form, ACQUIRE assumes 
that the resources of a VAX multi-tasking environment are available. 
For the convenience of users of more basic systems, ACQUIRE would need 
to be modified to include some simple processing capabilities. 

3.2 Simulations 

The effort that has gone into simulation has demonstrated the 
importance of planning, in the early stages of the program, the struc- 
ture and functions of the testing and validation methods to be used in 
the evaluation of the data-processing software. Since this program was 


initiated a great deal of knowledge and experience has been gained in 
these areas, and the need to devote adequate resources to the creation 
of a computer model reflecting as faithfully as possible the physical 
properties of the primary data is now clearly appreciated. Attempting 
to model the LTA system in all its material aspects would, however, have 
been an unacceptably labor-intensive task in the present program and 
some simplifications and approximations have been inevitable. For 
example, including all the properties of a full three-dimensional flow 
field would raise far too many complications, and the underlying assump- 
tion was therefore made that the flow is purely two-dimensional in the 
plane normal to the transmitter axis. Eventual extension to include 
three-dimensional effects on particle transits would, of course, be 
highly desirable. 



4.0 


RECOMMENDATIONS 


4.1 Software 

While the delivered software is complete and operational, our 
primary recommendation is to give priority to the program of testing, 
confirmation and verification in order to gain confidence that the 
programs really do what they are supposed to do, particularly with 

imperfect data. It remains necessary to inspect the correlograms 
visually to decide whether processing should continue or whether the 
data should be discarded. The criteria on which this choice is made, 
and which should be incorporated in the code have not been quantita- 
tively identified, and the degree of dependence on the experience of the 
operator remains higher than is really desirable. Progress towards 
reliable fully-automated software will depend on a thorough under- 
standing of the procedures in PROCESS and any shortcomings that may 

remain. Close technical liaison among users is also very important. 

In particular we recommend that PROCESS is exercised on high 
quality data over as large a range of values of the various parameters 
as possible, and only subsequently on progressively poorer data. The 
full test matrix is given in Table 4.1. 

The least satisfactory aspect of PROCESS in its present form is 
the underlying assumption that the flow can be modeled as a correlated 
bivariate Gaussian random process. This will certainly not be true in 
some cases - for example, in shear flows or where shock waves are 

present. A radical solution to this problem of general applicability 
would be to replace all the procedures following the initial condi- 
tioning of the correlograms (filtering, background removal and trunca- 
tion) with a model-free spline fit. 
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Table 4.1. Complete Test Matrix 


Spot Separation 
Spot Diameter 

852 micron 
40 micron 

Delay Time 
Number of Stores 

5 ns - 1 ys 
256 

Mean Velocity Components : 
u 

V 

2 - 1000 ms — 
2 - 1000 ms 

Turbulence Components : 
a 

X 

a 

v 

0 - 30% 
0 - 30% 

Correlation Coefficient 

-1 to +1 

Angular Steps 

3 

Angular Increment 

0.02° - 5° 

Height of Correlogram Peak, h^ 

10 - 104 

Fractional Background Level, h /h 

0.05 - 0.5 


n p 


A modification to the existing data-reduction procedure for Gaussian 
flows, which has rather more generality and which operates in (u,v) 
space, is recommended as a simpler solution and has the following 
s tructure: 


(i) 

(ii) 

(iii) 

(iv) 

(v) 

(vi) 

(vii) 


precondition correlograms in ( t, 9) space with present algorithms 
but in the order: filter, remove background, renormalize. 

transform to (u,v) space. 

find peak height H. at each angle 0 . and its coordinates 

fit Nth order polynomial (N = 4 may be a good choice) to the 

H. as a function of q. and find peak H and corresponding angle 

0^; this is taken to ile the mean flow Sirection. 
o 

fit Nth order polynomial to the set {(u^,Uj)} • 

find intersection of ray at 0 with this curve; this gives mean 
flow speed. 

find l//e points (or points at some other level) from correlo- 
grams in (u,v) space. 


(viii) perform elliptical fit as at present to find a , 

x 


a and p . 

y 

Implementing such a scheme would not involve major alterations to 


the structure of the present software. 



4.2 


Hardware 


Enhancements to the hardware which would improve the available 
accuracy in the estimation of turbulence could be made with existing 
technology. 

The first possibility is to modify the present Correlex by intro- 
ducing a known delay into the "A" channel. This has the effect of 
shifting the correlogram data to the left since the transit times are 
apparently smaller. The delay time per store may then be reduced to 
spread the data over a larger number of stores. Limitations on this 
technique are that the delay time cannot be less than 5 nanoseconds and 
that the "A 11 channel delay must be known to an accuracy consistent with 
that required for the measurements. 

Secondly, new correlators have been designed and constructed 
using a timing principle which is not based on discretization of time 
into set intervals. These instruments are capable of much improved 

accuracy and use a digital implementation of the former analog LTA 
filter/discriminators, together with improved algorithms. The pro- 
cessing of the event interval data is based on the Correlex principle. 
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